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Abstract Porcine deltacoronavirus (PDCoV) in Thailand 
was first detected in 2015. We performed a retrospective 
investigation of the presence of PDCoV in intestinal sam¬ 
ples collected from piglets with diarrhea in Thailand from 
2008 to 2015 using RT-PCR. PDCoV was found to be 
present as early as February 2013. Phylogenetic analysis 
demonstrated that all PDCoV variants from Thailand differ 
from those from other countries and belong to a novel 
group of PDCoV that is separate from the US and Chinese 
PDCoV variants. Evolutionary analysis suggested that the 
Thai PDCoV isolates probably diverged from a different 
ancestor from that of the Chinese and US PDCoV isolates 
and that this separation occurred after 1994. 

Porcine deltacoronavirus (PDCoV) causes an enteric dis¬ 
ease in pigs that is characterized by diarrhea and is clini¬ 
cally similar to porcine epidemic diarrhea (PED) and 
transmissible gastroenteritis (TGE) [4]. PDCoV, a recently 
emerged RNA virus of the family Coronaviridae , genus 
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Deltacoronavirus , was first reported in Hong Kong in 2012 
during an investigation to identify novel coronaviruses 
[14]. In February 2014, PDCoV was first detected in pigs 
with clinical diarrheal disease on farms in Ohio, USA. A 
subsequent retrospective investigation reported revealed 
that PDCoV had been present in the US as early as 2013 
[11]. PDCoV has since been identified in Canada, China, 
South Korea and Laos [1, 6, 9]. In Thailand, PDCoV was 
first reported in 2015 [10]. However, it is unknown whether 
PDCoV was present in Thailand prior to 2015. We there¬ 
fore conducted a retrospective study to detect the presence 
of PDCoV intestinal samples from Thailand that had been 
submitted for the identification of enteric diseases. Positive 
samples were subjected to full-length genome characteri¬ 
zation followed by genetic analyses to compare Thai 
strains of PDCoV to those from other countries and to 
investigate the genetic similarity and the evolutionary 
relationships among them. 

A total of 241 intestinal samples were assayed for the 
presence of PDCoV using RT-PCR. Intestinal samples 
were from 3- to 4-day-old piglets with clinical diarrheal 
disease collected from January 2008 to December 2015 and 
stored at -80 °C. The presence of other viruses causing 
enteric disease, including porcine epidemic diarrhea virus, 
and transmissible gastroenteritis virus was detected using 
RT-PCR. In brief, viral RNA was extracted from the 
samples using NucleoSpin RNA Virus (Macherey-Nagel 
Inc., PA, USA) in accordance with the manufacturer’s 
instructions. Complementary DNA (cDNA) was synthe¬ 
sized from the extracted RNA using M-MuLV reverse 
transcriptase (New England BioLabs Inc., Ipswich, MA, 
USA). PCR amplification was performed on cDNA using 
specific primers for the N gene of PDCoV, the spike gene 
of PEDV, and the N gene of TGEV as described previously 
[5, 10, 12]. 
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Samples positive for PDCoV were then subjected to full- 
length genome sequence characterization. Twenty-six pairs 
of primers were used to amplify the different regions of 
PDCoV [8]. PCR amplification was performed using Plat¬ 
inum® Taq DNA Polymerase High Fidelity (Invitrogen, CA, 
USA) in accordance with the manufacturer’s protocol. The 
PCR products were purified using a NucleoSpin Plasmid kit 
(Macherey-Nagel Inc., Bethlehem, PA, USA). The 5'- and 
3 / -terminal regions were determined by using a kit for rapid 
amplification of 5' and 3' cDNA ends (5' and 3'-RACE, 
Clontech, Japan). Sequencing was performed at First 
BASE Laboratory SdnBhd (Selangor, Malaysia) using an 
ABI Prism 3730XL DNA sequencer. 

The full-length PEDV genome sequence was assembled, 
and nucleotide and deduced amino acid sequences were 
aligned and analyzed using the CLUSTALW program [13]. 
The percentages of homology between the sequences at the 
nucleotide and amino acid levels were calculated. To 


investigate its evolutionary relationships to foreign PDCoV 
strains and estimate the divergence time of PDCoV in 
Thailand, phylogenetic analysis of the full-length nucleo¬ 
tide sequence of the PDCoV isolate collected in this study, 
together with 21 other PDCoV isolate sequences (Supple¬ 
mentary Table 1), was performed using the Bayesian 
Markov chain Monte Carlo (BMCMC) method imple¬ 
mented in the program BEAST vl.8.3 [2, 3]. A BEAST run 
was performed based on the GTR+G+I model with a 
coalescent Bayesian skyline tree prior and a clock model 
for each analysis using 200 million generations with sam¬ 
pling of every 10,000 generations and the first 10% dis¬ 
carded as burn-in. Tracer vl.6 was used to confirm that 
post-burn-in trees yielded an effective sample size (ESS) of 
>200 for all parameters. The resulting tree was viewed and 
generated in FigTree vl.4.2. 

The results of testing for PDCoV, PEDV, and TGEV in 
intestinal samples are shown in Table 1. Of 241 intestinal 


Table 1 Samples tested for porcine deltacoronavirus (PDCoV) by PCR in Thailand in 2008-2015 


Year 

Region 

Province 

Number of samples 

PDCoV positive (%) 

PEDV positive (%) 

TGEV positive (%) 

2015 

West 

Rachaburi 

24 

1 (4.17%) 

22 (91.67%) 

0 (0%) 



Nakhon Pathom 

14 

0 (0%) 

10 (71.43%) 

0 (0%) 


Middle 

Saraburi 

10 

1 (10%) 

10 (100%) 

0 (0%) 



Lopburi 

6 

2 (33.33%) 

6 (100%) 

0 (0%) 


East 

Chonburi 

11 

5 (45.45%) 

11 (100%) 

0 (0%) 



Chachoengsao 

2 

0 (0%) 

2 (100%) 

0 (0%) 


Northeast 

Nakhon Ratchasima 

8 

0 (0%) 

3 (37.5%) 

0 (0%) 



Buriram 

13 

0 (0%) 

13 (100%) 

0 (0%) 

2014 

West 

Ratchaburi 

12 

4 (33.33%) 

6 (50%) 

0 (0%) 



Nakhon Pathom 

3 

0 (0%) 

3 (100%) 

0 (0%) 


Middle 

Saraburi 

4 

0 (0%) 

4 (100%) 

0 (0%) 


East 

Chonburi 

1 

0 (0%) 

1 (100%) 

0 (0%) 

2013 

West 

Nakhon Pathom 

4 

0 (0%) 

4 (100%) 

0 (0%) 



Ratchaburi 

6 

2 (33.33%) 

4 (66.67%) 

0 (0%) 


East 

Chonburi 

7 

0 (0%) 

7 (100%) 

0 (0%) 

2012 

West 

Ratchaburi 

4 

0 (0%) 

4 (100%) 

0 (0%) 


Middle 

Saraburi 

2 

0 (0%) 

2 (100%) 

0 (0%) 


East 

Rayong 

1 

0 (0%) 

0 (0%) 

0 (0%) 

2011 

West 

Ratchaburi 

68 

0 (0%) 

9 (13.24%) 

0 (0%) 



Nakhon Pathom 

1 

0 (0%) 

0 (0%) 

0 (0%) 


Middle 

Saraburi 

9 

0 (0%) 

5 (55.56%) 

0 (0%) 


East 

Chonburi 

5 

0 (0%) 

1 (20%) 

0 (0%) 


Northeast 

Nakhon Ratchasima 

3 

0 (0%) 

1 (33.33%) 

0 (0%) 

2010 

West 

Ratchaburi 

9 

0 (0%) 

6 (66.67%) 

0 (0%) 


Northeast 

Nakhon Ratchasima 

2 

0 (0%) 

1 (50%) 

0 (0%) 



Udon Thani 

2 

0 (0%) 

1 (50%) 

0 (0%) 

2009 



1 

0 (0%) 

1 (100%) 

0 (0%) 

2008 

West 

Nakhon Pathom 

9 

0 (0%) 

7 (77.78%) 

0 (0%) 

Total 



241 

15 (6.22%) 

144 (59.75%) 

0 (0%) 
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Table 2 Gene lengths in four Thai PDCoV isolates and positions of 
compared to the other two PDCoV groups 

amino acid deletions/insertions and substitutions in 

Thai PDCoV isolates 

Thai isolates 


Gene 





ORFla/lb 

S 

E 

M 

N 

ST1_0213 Length 

Nucleotides 

18,786 

3,477 

249 

651 

1,026 


Amino acids 

6,262 

1,159 

83 

217 

342 

China group 

Deletions* 

401 LK 402 and 758 pvg 760 

571 v 

- 

- 

- 


Insertions* 

- 

51 n 

- 

- 

- 


Number of substitutions 

91 

35 s 

- 

1 

5 

US group 

Deletions* 

401 LK 402 and 758 PVG 760 

571 y 

- 

- 

- 


Insertions* 

- 

- 

- 

- 

- 


Number of substitutions 

96 

33 5 

- 

1 

5 

ST2_0214 Length 

Nucleotides 

18,786 

3,477 

249 

651 

1,026 


Amino acids 

6,262 

1,159 

83 

217 

342 

China group 

Deletions* 

401 LK 402 and 758 pvg 760 

571 v 

- 

- 

- 


Insertions* 

- 

51 n 

- 

- 

- 


Number of substitutions 

91 

37 5 

- 

1 

5 

US group 

Deletions* 

401 LK 402 and 758 pvg 7S0 

571 y 

- 

- 

- 


Insertions* 

- 

- 

- 

- 

- 


Number of substitutions 

96 

35 5 

- 

1 

5 

TT_1215 Length 

Nucleotides 

18,786 

3,477 

249 

651 

1,026 


Amino acids 

6,262 

1,159 

83 

217 

342 

China group 

Deletions* 

401 LK 402 and 758 pvg 760 

- 

- 

- 

- 


Insertions* 

- 

- 

- 

- 

- 


Number of substitutions 

90 

23 s 

- 

1 

4 

US group 

Deletions* 

401 LK 402 and 758 PVG 760 

- 

- 

- 

- 


Insertions* 

- 

- 

- 

- 

- 


Number of substitutions 

95 

23 5 

- 

1 

4 

NT1_1215 Length 

Nucleotides 

18,786 

3,479 

249 

651 

1,026 


Amino acids 

6,262 

1,160 

83 

217 

342 

China group 

Deletions* 

401 LK 402 and 758 pvg 760 

- 

- 

- 

- 


Insertions* 

- 

51 n 

- 

- 

- 


Number of substitutions 

139 

27 8 

- 

1 

3 

US group 

Deletions* 

401 LK 402 and 758 PVG 760 

- 

- 

- 

- 


Insertions* 

- 

- 

- 

- 

- 


Number of substitutions 

144 

27 8 

- 

1 

3 

* The numbers indicate the position 

in each gene 







Amino acid substitutions at different positions 


samples tested, all were negative for TGEV. Fifteen sam¬ 
ples (6.22%) were positive for PDCoV. One sample that 
was collected in February 2013 was PDCoV positive. The 
sample was from a farm in Ratchaburi, a province in the 
western region of Thailand where pig density is high. 
Interestingly, the sample was also positive for PEDV. Four 
samples collected in 2014 from two farms in Ratchaburi 
were PDCoV positive. On one of these farms, samples 
from 2013 and 2014 were positive for both PDCoV and 
PEDV. In 2015, PDCoV was detected in eight samples 
collected from five farms located in Ratchaburi, Saraburi, 


Fopburi and Chonburi. These four provinces contain three 
major swine-producing areas of Thailand. These results 
suggest that PDCoV is widespread throughout Thailand. It 
is noteworthy that 144 samples (59.75%) were positive for 
PEDV and that samples that were positive for PDCoV were 
also positive for PEDV. The results suggest the presence of 
mixed infections with PEDV and PDCoV in pigs. The 
production records of the herd in which PDCoV was first 
detected showed that there had been three diarrhea out¬ 
breaks during 2014 and 2015, and both PEDV and PDCoV 
were detected. In three diarrhea outbreaks, the pre-weaning 
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Table 3 Pairwise nucleotide and amino acid sequence identity values for the full-length genome sequence of the Thai PDCoV isolates from 
2013, 2014 and 2015 compared to PDCoV variants from China and the USA 

Year 

China PDCoV variants 


US PDCoV variants 



Nucleotide (%) 

Amino acid (%) 

Nucleotide (%) 

Amino acid (%) 

2013 and 2014 

97.2-97.8 

93.8-95.2 

97.3-97.4 

94.2-94.4 

2015 

96.9-97.8 

93.2-95.0 

97.0-97.3 

93.5-94.2 


mortality rate was 56.9, 52.1 and 65.3%, respectively. 
However, we were unable to define the relationship 
between PEDV and PDCoV regarding to the primary or 
secondary causes of the diarrhea outbreaks. Whether co- 
infection with these viruses would necessitate more-com¬ 
plex control methods or would result in more-frequent, 
repetitive outbreaks requires further investigation. 

Of the 15 PDCoV-positive samples, one from 2013, one 
from 2014, and two from 2015 were further characterized 


by full-length genome sequencing. The full-length genome 
sequences of the four Thai PDCoV isolates, one isolate 
detected in 2013 (P1_13_ST1_0213), one isolate detected 
in 2014 (P2_14_ST2_0214) and two isolates detected in 
2015 (P23_15_TT_1115 and P24_15_NT1_1215), are 
25,400-25,404 nucleotides (nt) in length, with a genome 
organization similar to that of all previously reported 
PDCoV genomes: 5’UTR-ORFla/lb-S-E-M-Nsp6-N- 
Nsp7-3’UTR [7, 14]. The ORFla/lb, S, E, M, Nsp6-N- 
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Fig. 1 Phylogenetic analysis of porcine deltacoronavirus (PDCoV) based on nucleotide sequences of the full-length genome performed using the 
Bayesian Markov chain Monte Carlo (BMCMC) method. The red dots represent Thai PDCoV isolates 
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Nsp7 genes are 18,786, 3,477, 249, 651, 282,1,026 and 600 
nucleotides in lengths, respectively (Table 2). Substitu¬ 
tions and deletions/insertions at the amino acid level 
occurred in several sequence regions, especially the 
ORFla/lb and S genes, of the four Thai PDCoV strains 
when compared to variants from the United States and 
China (Table 2). Interestingly, all four Thai PDCoV strains 
had amino acid deletions at 401 lk 402 and 758 pvg 760 . 

The pairwise nucleotide and amino acid sequence 
identity values of the full-length genome sequences of the 
four Thai PDCoV isolates are shown in Table 3. These 
sequences were 98.8-99.1% identical at the nucleotide and 
97.4-97.9% identical at the amino acid level and were more 
closely related to each other and to variants from China and 
the United States than to Thai strains isolated in 2015, 
suggesting that the Thai PDCoV population is continuously 
changing. 

Phylogenetic analysis based on the full-length genome 
sequence demonstrated that PDCoV isolates are divided into 
two groups, the Chinese- and US-like-groups, as was 
reported previously [1]. All four Thai PDCoV isolates were 
grouped in a novel cluster separate from variants from the 
United States and China but were more closely related to the 
Chinese isolates than to the US PDCoV isolates (Fig. 1). The 
currently circulating isolates of PDCoV share a common 
ancestor that was present around 1994, with 95 % HPD values 
covering 30 years from 1975 to 2004 (Fig. 1). Since then, 
Thai and US/China PDCoV clades diverged from the com¬ 
mon ancestor and evolved independently. The emergence of 
the Thai lineage was estimated to have occurred between 
2004 and 2012, while the US/China lineage evolved earlier, 
between 1994 and 2003. The clustering of the PDCoV iso¬ 
lates in the Thai PDCoV clade is well supported, with a 
posterior probability of 0.9993. Thai PDCoV sequences 
isolated in 2013 and 2014 grouped separately from Thai 
PDCoV sequences isolated in 2015. Moreover, all four Thai 
PDCoV isolates were descended from a common ancestor 
between 2004 and 2014 (Fig. 1). The analysis of amino acid 
changes also suggests that they are from a different lineage. 

In conclusion, we retrospectively investigated the presence 
of PDCoV in Thailand and found that PDCoV has been present 
in Thailand since at least February 2013 and that it is present in 
mixed infections with PEDV. Interestingly, the PDCoV vari¬ 
ants in Thailand are different from variants in other countries 
and belong to a novel group of PDCoV strains. Evolutionary 
analysis suggested that the Thai PDCoV isolates likely 
diverged from a different ancestor from that of the Chinese and 
US PDCoV isolates, and the separation occurred after 1994. 
This study provides interesting results associated with the 
emergence and the evolution of PDCoV. Further investigation, 
including more retrospective studies of available samples, will 
provide a better understanding of this pathogen. 
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